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Abstract 

We obtain a system for the spatial metric and extrinsic curvature of a spacelike 
slice that is hyperbolic non-strict in the sense of Leray and Ohya and is equivalent 
to the Einstein equations. Its characteristics are the light cone and the normal 
to the slice for any choice of lapse and shift functions, and it admits a well-posed 
causal Cauchy problem in a Gevrey class of index a = 2. The system becomes 
quasidiagonal hyperbolic if we posit a certain wave equation for the lapse function, 
and we can then relate the results to our previously obtained first order symmetric 
hyperbolic system for general relativity. 
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1 Introduction 



We consider, as in previous works 0, 0, the dynamics of General Relativity as the evo- 
lution of the geometry of spacelike slices of a spacetime manifold M x R. We take as 
coframe 

e° = dt, e* = dx l + pdt, (l) 

and write the spacetime metric 

ds 2 = -N 2 (6 f + g ij 9 i 9 j . (2) 

The geometry of a slice M t , x° = t, is determined by its metric g = g i jdx l dx : > and extrinsic 
curvature K^. The operator d , acting on time-dependent space tensors, is defined by 

d 

d o=gl~ C P> ( 3 ) 
where Lp is the Lie derivative with respect to the spatial shift vector (3 k . We have 

d o9ij = -2NK tj . (4) 

By using a combination fl^ = doRij — 2V(ji?j)o of the <% derivative of Rij with space 
derivatives of the momentum constraints, we have obtained a second-order system for 
Kij. This system reduces to a quasi-diagonal system for any choice of shift when the 
lapse iV is chosen so that either the time slicing is harmonic or the slices have a priori 
a given mean extrinsic curvature, H = K k k = h(x,t). The system obtained for (g, K) 
is hyperbolic and causal in the first case. It is a mixed hyperbolic (causal) and elliptic 
system in the second case, where one uses the i?oo equation to determine N. We have 
shown moreover that the system with the harmonic time slicing (and also generalizations 
of this slicing) can be written as a first order symmetric hyperbolic system if one uses 
also the i?oo equation. 

In this paper we use a combination of dod^R,-^ and derivatives of the constraints Riq 
and i?oo to obtain a system for (g, K) whose characteristics are the light cone and the 
normal to the spatial slice, for any choice of lapse and shift (as functions of spacetime). 
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This system is hyperbolic non-strict in the sense of Leray-Ohya[[3|. It admits a well-posed 
causal Cauchy problem in a Gevrey class of index a = 2. 

We also show that if the lapse N satisfies a wave equation whose source is an arbitrary 
(smooth) function of spacetime, g, K, N and first derivatives of N, the system becomes 
hyperbolic in the usual sense, with solutions in local Sobolev spaces. This system can 
be put in first order form, but it is not symmetric because the symmetrizing matrix is 
semi-definite. However, for an appropriate choice of the wave equation for the lapse, 
one can write our new system as a symmetric hyperbolic first order system that is the 
do derivative of the one previously obtained]!], [|. (We can also show that if the lapse 
is required to satisfy an appropriate elliptic equation, there is a well-posed hyperbolic- 
elliptic formulation M.) 



2 Non-strictly hyperbolic system 

We consider the combination 

Ay = d d Rij — 2<9 V(i-Rj) + VjVj-Roo = <9 f^ + ^V^/Zoo, (5) 

where (ij) = Mij + ji). We use the expression for fi^-, in the form given in [0] from the 
3+1 decomposition of the Ricci tensor, together with the expression for R 00 , to obtain 
the following identity 

Ay = doiNdK^ + doJij + ViVjiNdN-^K^K^+Cij, 

where 

□ = -N-%N-% + V fc V fc , (6) 

Ji, = d (HK, i:i -2K i k K jk ) + (N~ 2 d N + H)V l W j N 

-2iV- 1 (V fc iV)V (j (iV^ ) ) + 3(V fe A0V fc i^ (7) 
+N~ 1 K t;j V k (NV k N) - 2V {i (K j) k V k N) + N^HV^jN 2 
+2N- 1 (V {1 H)(V J) N 2 ) - 2NK k {l R j)k - 2NR kijm K km , 
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and 

dj = ViVjidoiN^doN) + Nd H) - BoiN-^iWjidoN + N 2 H)). (8) 

We see that contains terms of at most second order in K, and first order in g, after 
replacing dog%j by —2NKij. 

The identity given above shows that for a solution of the Einstein equations 

Raf3 = Pa/3, (9) 

the extrinsic curvature K satisfies, for any choice of lapse N and shift f3, a third order 
differential system which is quasidiagonal with principal part the hyperbolic operator 
<9oD. The other unknown g appears at second order except for third derivatives occuring 
in VjVjDiV. To avoid an explicit discussion of properties of matter evolution equations, 
we consider below the vacuum case. 

Theorem 1 Any solution of the vacuum Einstein equations R a/ 3 = satisfies the fol- 
lowing system of partial differential equations 

Ay = 0, d g i:i = -2NKij. (10) 

For any choice of shift [3 and lapse N > 0, (fJ7|) is a non-strictly hyperbolic system 
with domain of dependence determined by the light cone when the metric g is properly 
Riemannian. 

If the Cauchy data g(., 0) = 7, K(., 0) = k belong to the Gevrey class G^' 100 (Mo) , with 
7 properly Riemannian, there exists a solution of the system flT(\ ) in a neighborhood U of 
M , for any N, (3 G G2' loc (^); whose initial data <9oK(., 0) and 5o<9oK(., 0) are determined 
by the restrictions Rij(, .0) = and (doRij)(-,0) = 0. 

Proof. The principal matrix is triangular with diagonal elements either do or doO- 
Both of these operators are (strictly) hyperbolic, but the system is not (quasi- diagonal) 
hyperbolic, due to the presence of non-diagonal terms in the principal matrix. It is, how- 
ever, hyperbolic non-strict in the sense of Leray-Ohya . The Gevrey class in which the 
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system is well-posed is obtained by quasidiagonalization || and study of the multiplicity 
of the characteristics of the resulting diagonal elements. Here we obtain a quasidiagonal 
fourth order system for g, equivalent to the original system for g, K, by substituting 
Kij = — (2A^) _1 9 (7jj into A^. The principal operators on the diagonal are then d d n. 
They are hyperbolic non-strict with multiplicity two for Bq. The existence and uniqueness 
theorems of Leray-Ohya || give the result. 

3 Equivalence to the Einstein equations 

Theorem 2 Let N > and /3 be arbitrary in the Gevrey class G?, ,loc . Let (g, K) be a 
solution of the system 

Aij = Qij = dodopij - 2<9 V(i£j) + VjViPoo, (H) 

where p a p is some symmetric 2-tensor belonging to G^ Xoc that satisfies the conservation 
laws 

v a {p aP - \g aP g x »Px») = o. (12) 

Suppose the Cauchy data (g, K, d K, dod K.) Mo satisfy on M the Einstein equations 
(R a /3 — Pap) M = together with (do(Rij — pij)) m = 0. Then the metric (g, N, (3) satisfies 
the Einstein equations R a p — p a p = in the domain of dependence of Mo. 

Note. The equations (Rij — Pij)M = and (do(Rij — Pij))M — determine (doK^Mo 
and (B d Kij)M when the data (g, K) Mo are given. The equations (Ro a — Poo)m = 
are then equivalent to the usual constraints on (g, K)m , while the conditions (d (Ro a — 
Poo))m q = are a consequence of the restriction to M of the (twice-contracted) Bianchi 
identities and conservation laws. 

Proof. We set X af3 = R a/3 — p af3 . The Bianchi identities, together with the conservation 
laws, imply 

ViX'o - ^(N- 2 V X 00 + jtVaXij) = 0. (13) 
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We rewrite ( |TT|) 

VoVoXij - V V (i X,- )0 + VyViXoo = lhei, (14) 

where lhe n means linear and homogeneous in X and derivatives of order less than or 
equal to n. 

The sum of twice the Vo derivative of ( |1"3D and the contraction with of ([H]) gives 
the equation 

- iV- 2 V V X 00 + gVVjViXw = lh ei . (15) 
The Vq derivative of the Bianchi identity-conservation law with free index j reads 

V 2 V X ^ + VgViX« - V^V,(-iV- 2 Xoo + g mk X mk ) = 0. (16) 
We use (0) and the Ricci identity to obtain 

V 2 V X°^ + 2V V fc V (fc X^ - V j V k V k X 00 + iiV^V^Xoo (17) 

-Iv J [2V V fc X fe o- V k V k X 00 ] 
= Vo(-iV- 2 V 2 X^'o + V fc V A; X^)-^V J (-iV- 2 V 2 Xoo + V fc V fc Xoo) = lhe 2 . 
The equations (14), fll5D , (17), with unknowns X a p have a triangular principal matrix 



with diagonal principal operators □ and d^n. It is a non-strictly hyperbolic system. It is 
equivalent to the quasidiagonal system consisting of (jl~5|), flT6|) and the following equation 
obtained by combination of the previous ones 

V 2 (-iV- 2 V 2 + V k V k )X i3 = lhe 3 . (18) 

The system fllBD , (0), (|18[) is hyperbolic non-strict in the Gevrey class Gl' loc . The result 



follows from the Leray-Ohya uniqueness theorem. 

4 Hyperbolic system for (g, K, N) 

We add to the system (|TD|) a wave equation for N with arbitrary right hand side. Namely, 
we set 

DN = -N^doiN^doN) + V k V k N = F, (19) 



with F a given function of spacetime, N, g, K and first derivatives of N. Therefore, the 
equation = 0, reduced by the substitution of (|19|) and the replacement of d^g^ by 



—2NKij wherever it occurs, takes the form 

4n^i = /ii(2mg,K;3iniV). (20) 

The left hand side is third order in K and of order less than 3 in g and N. The numbers 
in fij denote the order of the highest derivatives occuring there. Because nN involves 
first derivatives of g, to obtain a hyperbolic system we differentiate (plj) by do and use 
&o9ij = —2NKij. This gives an equation which is third order in A and first order in g 
and K, 

B oN = d F. (21) 

Theorem 3 The system ([|) ; ft2(\), is a quasi- diagonal hyperbolic system. 

We attribute to the unknowns and equations the following Leray-Volevic indices: 

m(K) = 3, m(g) = 3, m(N) = 4, (22) 
n(20) = 0, n(|) = 2, n(0) = 1. 

The system assumes then a quasidiagonal form|| with principal operators 8q and dod. 
It is a hyperbolic system. 

We remark that a mixed hyperbolic-elliptic system for (g, K, N), with A^ determined 
by an elliptic equation is also possible ||. It is a straightforward extension of the result 
in M. 



5 First order system 

The system (^), (|T9|), ( |20|) can be put in first order form by a straightforward procedure. 
The system can be symmetrized, but the symmetrizing matrix is found to be semi- 
definite, not positive definite, so the usual existence and uniqueness theorems do not 
apply. It is an open question whether the system in first-order form is well-posed. 
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To achieve a symmetric hyperbolic system, the source F of the wave equation (|T9| ) 
for iV can be chosen so that no second spatial derivatives of K appear in the first order 
system: 

UN - NK mk K mk + NH 2 = F(g, N, b , a t ; x, t). (23) 

This equation is closely related to a form of harmonic time slicing that is allowed in our 
second order system for K , namely, 

d N + N 2 H = Nf(^a~\N;x,t), (24) 

where a(x, t) is an arbitrary positive scalar density of weight one and / is a scalar. Using 
fl2~4] ) and the Einstein equation for f?oo gives (^) with F = Fl^/ga^ 1 , N, 6 ; x, t), where 
bo = N^^^OqN is related to H by the harmonic condition fl24l) . In this case, the symmetric 
system obtained using ( |23"D is the time derivative of the first order symmetric hyperbolic 
system previously obtained in [T|. . 
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